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Abstract 

We present our first study of the time-dependent CP-violating asymmetry in B° — » J/tp7r° decays 
using e + e~ annihilation data collected with the BABAR detector at the T(4S) resonance during the 
years 1999-2002 at the PEP-II asymmetric-energy B Factory at SLAC. With about 88 million BB 
pairs, our preliminary results for the coefficients of the cosine and sine terms of the CP asymmetry 
are C J/V)7r o = 0.38 ± 0.41 (stat) ± 0.09 (syst) and Sj/^o = 0.05 ± 0.49 (stat) ± 0.16 (syst). 
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1 Introduction 



The Standard Model of electroweak interactions describes CP-violation in B meson decays by 
a complex phase in the three-generation Cabibbo-Kobayashi-Maskawa (CKM) [1] quark-mixing 
matrix. The b — > ccs modes that decay through charmonium, such as B° — ► J/t/iK®, yield precise 
measurements of the quantity sin2/3, where j3 = arg [ —V c ^* h /V t6 V£ h ] (see for example Refs. [2, 
3, 4]). The decay B° — > J/tpir is a Cabibbo-suppressed b —> ccd decay, whose tree contribution 
has the same weak phase as the 6 — > ccs modes (e.g. -B — ► J/ipK®). A portion of the penguin 
contribution has a different weak phase, which may give a time-dependent CP asymmetry that 
differs from the one observed in b — > ccs decays. 

In this measurement, about 88 million T(4S) — > -B-B decays are used to detect the decay chain 
5° — > J/tpir°, with J/0 — > e + e~ or J/t/j — > The BAR4R measurement of the i?° — > J/ipir 

branching fraction, (2.0 ± 0.6 (stat) ±0.2 (syst)) x 10 -5 , is described elsewhere [5]. Properties of 
the recoiling B meson are used to infer the flavor (B° or B°) of the B meson that is reconstructed 
from J/ip and 7r° candidates. The decay time distribution of B decays to a CP eigenstate with a 
B° or B° flavor ta g can be expressed in terms of a complex parameter A that depends on both the 
B°-B° oscillation amplitude and the amplitudes describing B° and B° decays to this final state [6] . 
The decay rate f+(f-) when the tagging meson is a B°(B°) is given by 



-\At\/r B() 



21m\ . . . 1 - |A| 2 

1 ± sm (Am d At) T cos (Am d Ai) 



(1) 



where At = t rec — t tag is the difference between the proper decay time of the reconstructed B 
meson (B rec ) and the proper decay time of the tagging B meson (Stag), T B° is the B° lifetime, 
and Amd is the B°-B° oscillation frequency. The sine term in Eq. 1 is due to the interference 
between direct decay and decay after flavor change, and the cosine term is due to the interference 
between two or more decay amplitudes with different weak and strong phases. Two amplitudes can 
contribute in the decay B° — > J/ip7r°. A portion of the penguin amplitude has the same weak phase 
as the tree amplitude, while the remainder of the penguin amplitude has a different weak phase. 
In b — > ccd decays, the tree contribution is Cabibbo-suppressed and the penguin and tree diagrams 
may enter at the same order, proportional to A 3 (where in this case A is the Wolfenstein parameter 
of the CKM matrix, rather than the complex parameter that appears in Eq. 1). Evidence for CP 
violation can be observed as a difference between the At distributions of B°- and /J°-tagged events 
or as an asymmetry with respect to At = for either flavor tag. We measure the two asymmetry 
coefficients, defined as 



2T rn\ , 1 - |A| 2 

where / is the final state. With these definitions, the absence of penguin contributions would 
give Sji^^o = — sin2/3 and Cj/^^o = 0. A statistically significant deviation from these values may 
indicate penguin contributions not only in B° — > J/V>7r°, but also in B° — > J/ipK® (at a reduced 
level governed by Cabibbo suppression). 
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2 The BABAR detector and dataset 



The data used in this measurement were collected with the BABAR detector at the PEP-II storage 
ring from 1999 to 2002. Approximately 81 fb -1 of e + e~ annihilation data taken at the T(4S) 
resonance are used, corresponding to a sample of about 88 million BB pairs. An additional 5 fb _1 
of data collected approximately 40MeV below the T(4S) resonance are used to characterize one of 
the background sources. 

The BABAR detector is described in detail elsewhere [7]. Surrounding the beam pipe is a 
silicon vertex tracker (SVT), which provides precise measurements of the trajectories of charged 
particles as they leave the e + e~ interaction point. A 40-layer drift chamber (DCH) surrounds 
the SVT, and both allow measurements of track momenta in a 1.5-T magnetic field as well as 
energy-loss measurements, which contribute to charged particle identification. Surrounding the 
DCH is a detector of internally reflected Cherenkov radiation (DIRC), which provides charged 
hadron identification. Outside of the DIRC is a CsI(Tl) electromagnetic calorimeter (EMC) that 
is used to detect photons, provide electron identification, and reconstruct neutral hadrons. The 
EMC is surrounded by the superconducting coil, which creates the magnetic field for momentum 
and charge measurements. Outside of the coil, the flux return yoke is instrumented with resistive 
plate chambers interspersed with iron (IFR) for the identification of muons and long-lived neutral 
hadrons. 

3 Candidate selection 

B° — > J/V>7T° candidates are selected by identifying J/ip — > e + e~ or J/ijj — > /t/ + /ii~ decays and 7r° 
— > 77 decays (details are given in Ref. [5]). For the J/ip — > e + e~ channel, photons consistent with 
bremsstrahlung are added and each lepton candidate must be consistent with the electron hypothe- 
sis. For the J/ip — > /i + channel, each lepton candidate must be consistent with the muon hypoth- 
esis. The invariant mass of the lepton pair is required to be between 2.95 and 3.14GeV/c 2 , and 3.06 
and 3.14GeV/c 2 , for the electron and muon channels, respectively. The photon candidates used to 
reconstruct the 7r° candidate are identified as clusters in the EMC within the polar angle range 
0.410 < #i a b < 2.409 rad that are spatially separated from every charged track, and have a mini- 
mum energy of 30MeV. The lateral energy distribution in the cluster is required to be consistent 
with a photon. The invariant mass of the photon pair is required to be 100 < m 77 < 160MeV/c 2 . 
Finally, the J/ip and tt° candidates defined above are combined using a mass-constrained kinematic 
vertexing algorithm. 

Two kinematic consistency requirements are applied to each B candidate. The difference, 
AE, between the B candidate energy and the beam energy in the center-of-mass frame must be 
—0.4 < AE < 0.4 GeV. The beam-energy-substituted mass, vtxes — ^/ (^^ cam ) 2 (.Pi?)^ i niust be 
5.2 < ttjes < 5.3 GeV/ c 2 , where -E^eam and p* B are the beam energy and B candidate momentum 
in the center-of-mass frame. 

Several kinematic and topological variables are linearly combined using a Fisher discriminant, 
to provide additional separation between signal and e + e~~ — > uu, dd, ss, cc (continuum) background 
events. The inputs to the Fisher discriminant are: the zeroth and second order Legendre polynomial 
momentum moments (Lq = Yi \Pi\ and L2 = Yli |Pi| 3cos 2 where p« are the momenta for the 
charged and neutral objects in the event that are not associated with the signal candidate, and 
9i are the angles between pj and the thrust axis of the signal candidate); the ratio of the second- 
order to zeroth-order Fox- Wolfram moment [8], computed using all charged and neutral objects not 
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Table 1: The efficiencies for the requirement on the Fisher discriminant and tagging, given inde- 
pendently, with statistical uncertainties. 



Source type 


Efficiency (%) of T > -0.8 


Tagging efficiency (%) 


B u -> J/ijjir u 


99.2 ±0.1 


65.6 ±0.6 


B° -> J/ipK^n ) background 


98.9 ±0.1 


65.6 ±0.6 


B — > J/ip X (inclusive J/ip) background 


94.9 ±0.7 


70.4 ± 1.4 


B — > X (SB generic) background 


98.5 ±0.4 


61.1 ± 1.6 


e + e~ — > (continuum) background 


28.6 ±0.7 


52.3 ±0.8 



associated with the signal candidate; | cos 0t | , where 0t is the angle between the thrust axis of the 
B candidate and the thrust axis of the remaining charged tracks and neutral objects in the event; 
| cos where 9i is the lepton helicity angle, defined as the angle between the negative lepton and 
B candidate directions in the J/ip rest frame. We require T > —0.8, which is 99% efficient for signal 
and rejects 71% of the continuum background. The efficiencies for satisfying this requirement are 
summarized in Table 1. 

4 Backgrounds 

The backgrounds come from decays which contain a J/ip particle or from purely random combina- 
tions. We split the backgrounds into four categories. 

One of the B — > J/ip X decays is B° — ► J/ip Kg(7r°TT°). In this case, one of the 7r°'s is emitted 
nearly at rest in the center-of-mass frame, and is thus missed in the reconstruction of the B 
candidate. The second is the more general class of B — > J/ip X (inclusive J/ip) decays, which 
contribute through random combinations of J/ip and 7T° candidates. This also includes cascade 
decays through other charmonium states, but excludes the specific B° — > J/tp Kg(7r°ir°) mode 
discussed above. Third is a purely combinatoric background contribution coming from the general 
decay B — > X (BB generic). Excluded from this definition are those decays already considered 
above. The fourth type of background is a combinatoric background due to u, d, s, and c quark 
production following the e + e~ annihilation, e + e _ — > qq (continuum). We study this background 
using an inverted lepton particle identification selection on the below-resonance data sample. In 
this case, the J/ip candidate is reconstructed from two particle candidates that are not consistent 
with a lepton hypothesis. Monte Carlo simulation is used to check that this procedure correctly 
models the background. 

5 Flavor tagging and measurement of At 

The methods for B flavor tagging, vertex reconstruction, and the determination of At, are described 
in Refs. [3, 9]. For flavor tagging, we exploit information from the recoil B decay in the event. The 
charges of energetic electrons and muons from semileptonic B decays, kaons, soft pions from D* + 
— ► D° 7r + decays, and high momentum particles are correlated with the flavor of the decaying B 
meson. For B decays, about 66% of the events can be assigned to one of four hierarchical, mutually 
exclusive tagging categories. The remaining untagged events are excluded from further analysis. 
The total tagging efficiency for each source type is shown in Table 1. 
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The time interval At between the two B decays is calculated from the measured separation Az 
between the decay vertex of the reconstructed B meson (B rec ) and the vertex of the flavor-tagging 
B meson (-Btag) along the beam axis (z axis). The calculation of At includes an event-by-event 
correction for the direction of the B rec with respect to the z direction in the T(4S) frame. We 
determine the z position of the B rec vertex from the reconstructed vertex of the J/ip candidate. 
The S tag vertex is determined by fitting the tracks not belonging to the B vec candidate to a common 
vertex. An additional constraint on the tagging vertex comes from a pseudo-track computed from 
the -B rec vertex and three-momentum, the beam-spot (with a vertical size of 10 /im), and the T(4S) 
momentum. For 99.5% of the reconstructed events the r.m.s. Az resolution is 180 fim. Convergence 
is required for both the B rec and B tag vertex fits. Finally, At must be between —20 and 20 ps, and 
it is required to have an uncertainty satisfying <7 At < 2.4 ps. 



6 Maximum likelihood fitting technique 

We extract the CP asymmetry by performing an unbinned extended likelihood fit. The likelihood 
is constructed from the probability density functions for the discriminating variables iues, AE 1 , 
and At. The quantity that is maximized is 



e-^U^ N 

C = 



nff Si 9 „„. -p Sig V Sig T> Sig 
Ua n Sig I toes AE 'At 

i fKs T)Ks T>Ks 

+ Jot n Ks /-TOES-A.E ' At 

I fine q-)Inc q-jlnc 

~rj a "Inc ' m^-AE 'At 

+ fa B n BB ^TOES ^a! ^At 
, fCont _ T>Cont nC'ont >j-)Cont~\ 

+ J a n C ont /''toes ' AE < J At h 



(3) 



where rij is the number of events for each of the 5 hypotheses (1 signal and 4 backgrounds) and N is 
the number of input events. The V are the probability density functions (PDFs) for each discrim- 
inating variable and signal or background type. The parameters fa are the tagging efficiencies for 
each of the 4 tagging categories a and each of the signal or background types j. For the B° — > J/ip n° 
signal and B° — > J/tp K®(Tr ir°) background, the values of f a are measured with a sample (-Bfl av ) of 
neutral B decays to flavor eigenstates consisting of the channels D^*^ h + (h + = ir + , p + , and af) and 
J/ijK*°(K*° -> K+-K-) [3]. For the inclusive J/tp background and BB generic background, they 
are measured with Monte Carlo simulation [10], and for the continuum background, they are mea- 
sured with the inverted lepton particle identification data sample. We discuss the discriminating 
variables (toes, and At) in the following sections. 

6.1 Probability density functions for m ES and AE 

The signal toes distribution is modeled as the sum of two components. The first is a modi- 
fied Gaussian function which, for values less than the mean, has a width parameter that scales 
linearly with the distance from the mean. The second component, accounting for less than 6% 
of the distribution, is an ARGUS function [11], which is a phase-space distribution of the form 
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m ES\/(l ~~ ^ ES ) exp(£(l — ^ ES )), with a kinematic cut-off at Eb eam = 5.289 GeV, and one pa- 

V beam beam 

rameter to fit in the exponential, £. The signal AE distribution is modeled by the sum of a Crystal 
Ball function [12] and a second order polynomial. The Crystal Ball function is defined as 



C(AE) = < 



e 2^ if AE > m — aa 

( — ) n e~~2~ 

^ r -^i; ar if AE < m — aa, 



(4) 

where m is the position of the maximum. The parameter a determines the cross-over point from a 
Gaussian behavior to a power-law, and is expressed in units of the peak width a. The parameter 
n is a real number which enters into the power-law portion of the function. The parameters of 
these PDFs are determined by fitting to a signal Monte Carlo sample. The maximum of the 
AE distribution is a free parameter of the full CP likelihood fit to allow for EMC energy scale 
uncertainties. 

The kinematic variables 7tt.es and AE are correlated for the B° — > J/ip and the 

inclusive J/ip backgrounds. To account for this, 2-dimensional PDFs are employed. Variably 
binned interpolated 2-dimensional histograms of these variables are constructed from the relevant 
Monte Carlo samples. 

The 777.es PDFs for the BB generic and continuum backgrounds are modeled by the ARGUS 
function and the AE PDFs for these two backgrounds are modeled by a second order polynomial. 
The parameters for these PDFs are obtained from the BB generic Monte Carlo sample and the 
inverted lepton particle identification data sample. 

6.2 Probability density functions for At 

The PDFs used to describe the At distributions of the signal and background sources are each a 
convolution of a resolution function 7Z and decay time distribution V: V(At) = TZ(At) <8) V(At). 

For the signal, the resolution function [9] consists of the sum of three Gaussians, which will be 
referred to as the core, tail, and outlier. The means of the Gaussians are biased away from zero 
due to the charm content of the side of the event used for tagging. For the core and tail Gaussians 
this bias is multiplied by the At per-event error a At- The widths of the core and tail Gaussians are 
the products of the per-event errors and scale factors. The tail Gaussian has a fixed scale factor of 
3 and the outlier Gaussian has a fixed width of 8 ps and zero mean. The five remaining parameters 
are measured with the large -Bfl av data sample. The bias of the core Gaussian has different values 
for each of the four tagging categories. 

The decay time distribution is given by Eq. 1 modified for the effects of B flavor tagging: 

e -|At|/T fl0 

Z£ f (At) = — {(1 T Aw a ) ±S f (l- 2w a ) sm(Am d At) 

" 4r go 

T C f (1 - 2w a ) cos(A777 d At)}, (5) 

where T>+ f CD~ f ) is for a B°(B°) tagging meson. The variable w a is the average probability of 
incorrectly tagging a B° as a -B (wf ) or a B° as a B° {w^°), and Aw a = — . Both w a 
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Table 2: Results of the CP likelihood fit. Errors are statistical only. The global correlation is 
0.14 for Cj/^^o and 0.15 for Sj^ w o. The projections of the PDFs are shown in Figure 1 and the 
asymmetry in Figure 2. 





Fit results 


Cj/l/>7T° 


0.38 ±0.41 




0.05 ± 0.49 


Signal AE Maximum ( MeV) 


-13.2 ±7.2 


i? u — ► J/ij)^ signal (Events) 


40 ±7 


B° -> J/t/jK^i^TT ) background (Events) 


140 ± 19 


Inclusive J/ip background (Events) 


109 ± 35 


BB generic background (Events) 


52 ±25 


Continuum background (Events) 


97 ±22 



and Aw a are determined using the i?fj av data sample [3] . The values of Am^ and t b o are the 2002 
PDG averages [13]. 

The PDF used to model the At distribution for the B° — > J/ip Kg(ir°7r°) background takes the 
same form as that for signal, but with Sj^ K o = 0.75 [14], and Cjj^ K o = 0. 

The parameterizations of the At PDFs for the inclusive J/ip background and the BB generic 
background consist of lifetime and prompt components. The resolution function for each component 
is the sum of core and outlier Gaussians, where the width of each core Gaussian is the product of 
a At and a scale factor. Once again, the width and mean of the outlier Gaussians are fixed to 8ps 
and zero respectively. For each of these background sources, the fraction which is in the lifetime 
component, the decay lifetime parameter, and the resolution parameters are the values determined 
from the Monte Carlo simulation. 

The At PDF for the continuum background consists of a double Gaussian which has the same 
form as the prompt component of the inclusive J/ip and BB generic At PDFs, where in this case 
the parameter values are obtained by fitting the inverted lepton particle identification data sample. 

6.3 Results of the CP asymmetry fit 

The results of the CP asymmetry fit to 438 events found in 81 fb _1 of data are shown in Table 2. 
The projections in toes, AE 1 , and At, are shown in Figure 1. Figure 2 shows the At distributions 
and asymmetries in yields between B and B° flavor tags as functions of At, overlaid with the 
projection of the likelihood fit results. 

7 Systematic uncertainties 

The contributions to the systematic errors in Cj/^^o and Sj/^^o are summarized in Table 3. The first 
class of uncertainties are those obtained by variation of the parameters used in the toes, AE, and At 
PDFs, where the dominant sources are the uncertainties in the signal AE PDF parameters. Another 
contribution is due to the energy scale uncertainties in the modeling of the B° — ► J/ip K®(ir°ir°) 
background. An additional systematic uncertainty comes from altering the configuration of the 2- 
dimensional PDFs for the B° -> J/^K°(7r 7r° ) and inclusive J/ip backgrounds. A systematic error 
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Figure 1: Projections in a) mEs> 

c) AE, and e) At for the results of the CP fit to 81 ftr 1 of 
data. The legend in a) applies to the plots on the left hand side. The projection in b) uies is 
shown with the requirement —0.11 < AE < 0.11 GeV. The projection in d) AE is shown with 
the requirement thes > 5.27 GeV/ c 2 . The projection in f) At is shown with the requirements 
—0.11 < AE < 0.11 GeV and m^s > 5.27 GeV/ c 2 . These plots do not represent the full information 
used in the maximum likelihood fit, but only a partial view of the data. 
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A t (ps) 




Figure 2: Number of candidates in the signal region a) with a B° tag N B o and b) with a B° tag 
iV^o, and c) the raw asymmetry (N B o — N^ )/(N B o + N B o), as functions of At. Candidates in 
these plots are required to satisfy —0.11 < AE < 0.11 GeV and tues > 5.27GeV/c 2 . The curves in 
a) and b) are projections that use the values of the other variables in the likelihood to determine 
the contribution to the signal or one of the backgrounds. 
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Table 3: Summary of the systematic errors. 



Source 


Error on Cj/^o 


Error on Sj/^o 


Parameter variations 


ttjes an d AE" parameters 


0.048 


0.130 


Tagging fractions 


0.002 


0.007 


At parameters 


0.027 


0.022 


Additional systematics 


EMC energy scale B u -» J/^^(vr u 7r u ) 


0.009 


0.002 


Changing the 2-D histogram PDFs 


0.009 


0.029 


Using 2-D PDF for signal 


0.073 


0.079 


Beam spot, boost/vtx., SVT misalign. 


0.012 


0.012 


Total systematic uncertainty 


0.093 


0.157 



to account for a correlation between the tails of the signal ttjes and AE distributions is obtained 
by using a 2-dimensional PDF. 

8 Summary 

An unbinned extended maximum likelihood fit has been performed on 81 fb _1 of data collected at 
BABAR, yielding preliminary values for the coefficients of the cosine and sine terms of the time- 
dependent CP asymmetry in B° — ► J/tpir decays: 

Cj/^o = 0.38 ± 0.41 (stat) ± 0.09 (syst), 
SjtyTfl = 0.05 ± 0.49 (stat) ± 0.16 (syst). 
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